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Abstract 

We study the T exclusive decay into double charmonium, specifically, the 5-wave charmonium 
J/tp plus the P-wave charmonium XcO l 2 in the NRQCD factorization framework. Three distinct 
decay mechanisms, i.e., the strong, electromagnetic and radiative decay channels are included 
and their interference effects are investigated. The decay processes T(15, 2S, 3S) — > J/tp + Xci,o 
are predicted to have the branching fractions of order 10 -6 , which should be observed in the 
prospective Super B factory. 
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I. INTRODUCTION 



Anyone who has ever browsed the Meson Summary Table in the biannual review of 
particle data group will be impressed by the extremely rich decay channels of the B, D, 
and J/ip mesons Unlike the heavy flavored mesons, which can only decay via the 
weak interaction, the unflavored heavy quarkonia, decay through the heavy quark- ant iquark 
annihilation initiated exclusively by the strong and electromagnetic interactions. Although 
quite a few decay channels have been established for the charmonia system over the past 
few decades, the experimental information about the bottomonia decay is still very sparse. 
Until very recently, the Belle Collaboration has observed a few exclusive decay channels 
of T(1S, 2S) into light hadrons for the first time, e.g., into the vector-tensor states and the 
axial- vector-pseudoscalar states j2[. 

Because of the much more copious phase space opened at the bottomonium energy level, 
the typical branching fraction for a given hadronic decay mode of a bottomonium is greatly 
diluted with respect to that of a charmonium. Thanks to the weaker strong coupling at the 
bottom mass scale, the perturbative QCD is expected to work more reliably for the hadronic 
bottomonium decay than for the charmonium. 

An interesting class of hadronic decay processes of bottomonium is into double charmo- 
nium, which may presumably be predicted with less uncertainties than into the light hadrons. 
In the recent years, some exclusive decay processes of bottomonium into double charmonia 
have been intensively studied in the perturbative QCD framework, e.g., i] b — > J/ip J/ip 



Xw,i,2 J/^J/ip |7Hll|. and T — > J/ip + rj c [12|, |l3j. These studies are largely inspired 
by the various double-charmonium production processes in e + e~ annihilation, which were 
first observed at the B factories a decade ago |l4j-|l6|. Triggered by the disquieting discrep- 
ancy between data and theory, since then a great number of theoretical studies have been 
conducted for the processes e + e~ — > J/ip + i] c 17-29[, and e + e~ — >■ J/ip + Xc0,i,2 (30-32] . 



Besides the tremendous amount of data near the T(4S), the Belle experiment to date 
has also collected about 102 million T(IS') samples and 158 million T(2S) samples. There- 
fore, it appears more promising to observe the double-charmonium production from the T 



decay than from the C-even bottomonia decay. In Ref. [13(|, the exclusive decay of T to 



a vector-pseudoscalar charmonium states was studied in the framework of nonrelativistic 



QCD (NRQCD) factorization 33]. The corresponding branching fraction was estimated to 
be of order 10~ 6 and seems to have a good chance to be observed at the Super B factory. 
In this work, we further investigate the T decay into the J/ip plus a spin-triplet P-wave 
charmonium Xcj (.7=0,1,2). This work should be considered as a sequel of Ref. [l3| 1 . 

Although neither of these exclusive decay modes have been observed yet, several upper 
bounds for T(nS) inclusive decay into J/ip or Xcj have already been placed experimen- 
tally (J: 

B[T(1S) ->• J/ip + X} = (6.5 ± 0.7) x 10" 4 , B[?(1S) X co + X] < 5 x 10~ 3 , 
B[T(1S) -»• Xci +X] = (2.3 ± 0.7) x 10~ 4 , B[T(1S) ->• Xc2 + X] = (3.4 ± 1.0) x 10" 4 
B[T(2S) -)• J/ip + X] < 6 x 10~ 3 , B[T{AS) -> J/ip + X] < 1.9 x 10" 4 . (1) 



1 The main results of this paper have already been reported in Ref. 3J|. Nevertheless, some significant 



im pro vements have been made in the current work: i.e., some errors in calculating the three-gluon channel 



34j have been corrected, and the contribution from the two-gluon-one-photon channel is also included. 
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It will be interesting to examine to which extent these upper bounds are saturated by the 
predicted branching fractions for T(nS) — > J /ip + Xc0,i,2- 

As is well known, the hadronic decay of T can be categorized into three distinct classes: 
bb either annihilating into three gluons (strong decay), or single photon (electromagnetic 
decay), or two-gluons and a photon (radiative decay). On the experimental ground, the 
inclusive decay rates from these three decay channels have been available long ago [l[: 

B[T -> ggg\ : B[T -> 7* -)• X] : B[T ggj] = 82.7% : 7.5% : 2.2%, (2) 

where these three branching ratios sum up to 1 — B[T — > = 92.5%, as they should 2 . 

For the exclusive hadronic decay T — > J/tj) + Xco,i,2, one may also be interested in as- 
certaining the relative strength and the interference pattern among these different decay 
channels. This sort of study has been conducted for the process T — > J ftp + rj c [l3| . Note 
that there have lasted constant experimental efforts to infer the relative phase between 
the strong and electromagnetic amplitudes in exclusive J/ip and ip' decays into two light 
mesons |35H40J. As we will see later, being a consequence of nib > m c ^> Aqcd, the relative 
phases among three distinct decay channels in our processes stem from the short-distance 
loop contribution, which can actually be calculated in perturbation theory. 

The rest of the paper is organized as follows. In Section [TTJ we express the polarized 
and unpolarized decay rates in terms of the helicity amplitudes and briefly state the helicity 
selection rule. In Section [TTT| we conduct the lowest order (LO) calculation for each indepen- 
dent helicity amplitude associated with the decays T — > J/ijj + Xco,i,2, within the NRQCD 
factorization approach. The contributions from three distinct decay channels, i.e., electro- 
magnetic, strong, and radiative decay channels, are all included, and the analytic expressions 
for each of the helicity amplitudes are given. In Section HV] we present our predictions of the 
interference pattern among three distinct decay channels for T — > J/tp + Xco,i,2, and of the 
polarized and unpolarized partial decay widths and the corresponding branching fractions 
for T(1S,2S,3S) decays into J/ip + Xco,i,2- We find that it appears quite promising for 
the prospective Super B experiment to observe these hadronic decay processes. Finally we 
summarize in Section |V] In Appendix |AJ we list the explicit expressions of the 10 helicity 
projectors that are used in Section IIII1 



II. POLARIZED DECAY RATES AND HELICITY SELECTION RULE 



It is of some advantage to utilize the helicity amplitude formalism [41l . |42| to analyze the 
hard exclusive reactions, in particular for the decay process studied in this work. From the 
experimental perspective, the helicity amplitudes can in principle be accessed via measuring 
the angular distributions of the decay products of J/tp and Xcj (J = 0,1,2), provided 
that the statistics is sufficient. From the theoretical viewpoint, some essential dynamics 
underlying perturbative QCD is clearly encoded in the helicity amplitude analysis, which 
becomes rather obscured if one only looks at the unpolarized reaction rates. 

We will work in the T rest frame throughout this work. Suppose the spin projection of 
the T along the z axis to be S z (The z axis, say, may be chosen as the beam direction of 
the e~ and e + collider which resonantly produce a T meson). Let A, A denote the helicities 



2 We have not included the contribution from the radiative transition T — > 77^7, which has a completely 
negligible branching ratio. 
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carried by the outgoing J /if) and Xcj, respectively, and 6 signify the angle between the 
J /ib momentum P and the z axis. The differential polarized decay rate can be expressed 
as gj 52 



dT[?(S z )^J/iP(\) + X cj(\)} 
dcos8 



IQtiM 2 



4,a-#) 



A 



A. A 



(3) 



where A J X r (J = 0, 1,2) characterizes the corresponding helicity amplitude, which encom- 



passes all the nontrivial QCD dynamics. The angular distribution is fully dictated by the 
quantum numbers S z , A and A through the Wigner rotation matrix d 3 mm ,(9). Note that 

angular momentum conservation constrains that |A — A| < 1. In ([3]), the magnitude of the 
three- momentum carried by the J/ib (or Xcj) is determined by 



AV 2 (M|,M^,M x 2 ej 
2M T 



(4) 



where A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. 

Integrating ([3]) over the polar angle, and averaging over all three possible T polarizations, 
one finds the integrated rate of T decay into J/ib + Xcj in the helicity configuration (A, A) 
to be 



r[T ->■ J/ib{\) + Xcj(A)] 
|P| A'- 

A, A 



16%M 2 



A 



A. A 



J_d cos 9 ^J2\ d k,x-x 



(0) 



2\-kM\ 



(5) 



Since this decay process can be initiated by the strong or electromagnetic interactions, one 
can resort to the parity invariance to reduce the number of independent helicity amplitudes: 



A 



J A J 



A. A 



■lYA 



-A,-A" 



(6) 



As a consequence, the helicity channel T — > J/ib(0) + x c i(0) is strictly forbidden. 

Starting from ([5]), one readily obtains the unpolarized decay rate by summing the con- 
tributions from all the allowed helicity channels: 



r[T^ J/ip + Xd] 
r[T^ J/ip + Xc2] 



24nM 2 



2AnM 2 



24nM 2 



A I 

^0,0 1 



2 \A° I 

1*^-1,0 | 

^2 |*/1.-|^q| -)- 2 |^4g^| -)- 2 \A--y ]J 



(7a) 
(7b) 



2 |2 
0,0 1 



2 ^ * + 2 1*4^ + 2 \A{^ + 2 \A{ >2 \* (.7c) 



There are two, three and five independent helicity amplitudes for T — > J/ib + Xcj (J = 
0,1,2), respectively, as enforced by the angular momentum conservation. We have also 
included a factor of 2 to account for the parity-doublet contributions. 

One important piece of physics underlying the hard exclusive reactions is that each helicity 
amplitude possesses a definite power-law scaling in the inverse power of large momentum 
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transfer, controlled by the helicity selection rule (HSR) [43]. At asymptotically large m&, 
the polarized decay rate in our process scales as [171 ]: 

r[T ^ j/^(A) + XcJ (A)] 

here t> denotes the characteristic velocity of the charm quark inside a charmonium. 

Equation ([H]) implies that the helicity state which possesses the slowest asymptotic de- 
crease, is the one that conserves the hadron helicities |A + A] = 0. In line with the angular 
momentum conservation, the only possible configuration is (A, A) = (0,0). For each unit of 
the violation of the helicity conservation, there is a further suppression factor of l/mf . In 
the limit nib oo, perhaps only the (0, 0) helicity state is phenomenological relevant. Note 
that in NRQCD factorization language, the charm quark is also treated as heavy, and in 
fact its mass acts as the agent of violating the hadron helicity conservation. 

We note that, the power-law scaling specified in ([8]) is in general subject to the mild 
modifications due to the ln(m^/m^) from the loop contribution. This logarithmic scaling 
violation will be examined in detail in Section IIHI 




III. THE CALCULATION OF THE HELICITY AMPLITUDES IN NRQCD FAC- 
TORIZATION APPROACH 



The hard exclusive decay process T — > J/ip + Xcj (J = 0, 1,2) is characterized by two 
hard scales set by the bottom and charm quark masses. This process can proceed via three 
separate channels: the bb pair first annihilates into a single photon, or three gluons, or two 
gluons plus a photon, subsequently the highly virtual photon/gluons transition into two cc 
pairs, which finally materialize into two fast-moving charmonium states. 

Two influential perturbative QCD approaches are legitimate to describe such type of 



decay process, i.e., the light-cone approach [441 . 145| which is based on twist expansion, and 
the NRQCD factorization approach [33] that is based on the quark velocity expansion. As 
was seen in Sec. HIl most helicity channels associated with the process T — > J/ip + Xco,i,2 are 
of helicity-suppressed type. This feature impairs the practical usefulness of the light-cone 
approach, since the higher-twist light-cone distribution amplitudes of charmonia are rather 
poorly understood at present. 

On the other hand, the NRQCD factorization approach, which is based upon a completely 
different expansion strategy, does not confront any obstacle in dealing with helicity-flipped 
channels. In the past two decades, this framework has been widely applied to numerous 
quarkonium decay and production processes jsjjj]. In contrast with the light-cone approach, 
the nonpertubative input parameters in NRQCD factorization approach are numbers (local 
NRQCD matrix elements, or wave functions at the origin) rather than functions (light-cone 
distribution amplitudes). In this regard, NRQCD approach seems to be more economic and 
predictive than the light-cone approach. 

In this work, we will investigate the process T — > J/ip + Xco,i,2 i n the framework of 
NRQCD factorization, incorporating aforementioned three distinct decay mechanisms 3 . We 



Unlike the exclusive double-charmonium production in e + e~ annihilation, where a factorization theorem 
in NRQCD has been proved to all orders in a s 461 ] ■ there has not yet existed any rigorous proof for the 
validity of NRQCD approach to T — > J/tjj + Xc0,i,2- 
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will be content with the lowest order accuracy in both the velocity expansion and the strong- 
coupling constant expansion. We are aware that our results may be subject to considerable 
uncertainty from various sources, yet still hope the predicted decay rates may capture the 
correct order of magnitude. 

At the LO in the bottom and charm velocity, one can expedite the NRQCD approach 
calculation by invoking the covariant projection method [l7|, i.e., first calculate the on- 
shell T-matrix for bb(Q) — > cc(P) + cc(P), then project each quark-antiquark pair onto the 
intended spin-color and orbital angular momentum states. In this case, all the involved 
nonpertubative quarkonium-to-vacuum NRQCD matrix elements can be well approximated 
by three (the first derivative of ) wave functions at the origin for the quarkonia T, J/i/j, 
and Xcj'- -Rx(O), -Rj/^,(0), and R' Xcj (0). Each of them can be either obtained from the 
quark potential models, or calculated from lattice simulation, or directly extracted from the 
quarkonia decay data. 

The product of these three nonpertubative wave functions at the origin ubiquitously 
enters each helicity amplitude. Thus it seems convenient to define a reduced dimensionless 
helicity amplitude, of which these nonperturbative factors are explicitly pulled out. First, 
let us introduce a mass ratio variable: 



m 



r = ( 9 ) 



The reduced helicity amplitude, dubbed a^~, is related to the standard helicity amplitude 



as follows: 



. — 3 R r )R JN) )R' ^n^t, , 
' JW h mi > v ; 



where N c = 3 denotes the number of colors. Note that the scaling factor dictated by the 
HSR has been explicitly factored out, the reduced amplitude a^- is thereby expected to 

scale with r as O(r ). 

Inserting (jTUl) back into ([5]), one can reexpress the integrated polarized decay rate as 

IPI 



r[T jma) + XcJ (X)} = iv c 3 ^(o)^(o)^ 2 eJ (o) 3M 2 rr ; bm ^ 2+|A+A| 



x 



(11) 



for each helicity channel. The subscripts 7, 3g, and 'ygg emphasize the decay channel with 
which the reduce amplitude is affiliated. Obviously, it is of interest to ascertain the relative 
strength and phase among these different types of amplitudes. 

In the remainder of this section, we will present the analytic expressions of the reduced 
helicity amplitudes associated with each decay channel. 



A. Single-photon channel 



We start by considering the decay channel T — > 7* — > J ftp + Xcj, with some typical 
LO diagrams shown in Fig. [lj This process is very similar to the continuum J ftp + Xcj 



production in e + e annihilation [17 
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FIG. 1: Two representative lowest-order diagrams that contribute to T — > 7* — > J/ip + Xcj- There 
are totally four diagrams in class a) and two in class b). 



After obtaining the decay amplitude AZ in NRQCD factorization, one can employ the 
helicity projectors enumerated in Appendix |A]to project out 10 corresponding helicity ampli- 
tudes. It is straightforward to follow equation ( ITU]) to read off the reduced helicity amplitude 
in the single-photon channel: 



7;A,A 



-e b e c aa s c'l-Ar) 



(12) 



where e\, — — | and e c = | are the electric charges of the b and c quarks, a and a s are the 
QED and QCD coupling constants, respectively. The coefficient functions c^(r) read 



-0,0 v 
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c o,o( r 
c lo( r 



1 + 10r - 12r 2 + 2ry 



9-14r--(--6 
2 \r 



'13a) 



-^6 
■2VE 



2 - 7r + -y 
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y/2[l - 2r - 12r 2 + 2ry] 
V2 



V6 



2\r 
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(13b) 



V3 



3 - llr - - ( - - 3 
2 \r 



2- V - 

r J 



c? x (r) =2V6 



y ( 1 



4 V r 



(13c) 



where y = —a/a s . The y-dependent terms characterize the photon fragmentation contribu- 
tions as depicted in Fig. Wp), which are often accompanied by an enhancement factor 1/r 
for the transversely-polarized J/ip. 

Barring the pure QED fragmentation contributions, these 10 coefficient functions agree, 
up to an immaterial phase, with those associated with the process e + e _ — > J/ip+Xco,i,z 4 - 



4 We take this opportunity to point out a typo in equation (8c) in Ref. [32[, where (r) was erroneously 
typed as y/2(ll — 3r). 
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FIG. 2: Some representative lowest-order diagrams that contribute to T — >• 3g — > J/ip Xcj- 



It is interesting to mention that, for some accidental reason, the QCD part of the single- 
photon J/ip(±l) + Xci(O) amplitude receives an extra suppression factor than implied from 
HSR. 



B. Three-gluon channel 

Next we turn to the strong decay channel T — > 3g — > J/ip + Xcj, which supposedly 
makes the most significant contribution. Some of the representative LO diagrams have been 
illustrated in Fig. [2j In contrast with the single-photon channel, this channel first starts at 
the one-loop level. The charge conjugation invariance guarantees that one needs only retain 
those diagrams with "Abelian" gluon topology as shown in Fig. [2j 

After obtaining the decay amplitude A^ g using the covariant projection technique, and 
prior to performing the loop integration, we apply those helicity projection operators given 
in Appendix |A] to project out 10 corresponding helicity amplitudes. This operation brings 
forth great simplification, because all the polarization vectors (tensors) of T, J/ip and Xcj 
have been eliminated from the integrand, and the numerators in loop integrals now become 
Lorentz scalars comprised entirely of the external and loop momenta. 

It is then straightforward to utilize the partial fractioning technique to reduce all the 
higher-point one-loop integrals into a set of 2-point and 3-point scalar integrals. Most of the 
encountered scalar integrals can be found in the appendix of Ref. [13|, whose correctness has 



been numerically checked by the Mathematica package LoopTools [42]. There also arise 
some nonstandard 2- and 3-point scalar integrals, which contain propagators with quadratic 
power due to the projection of the P-wave state. All of their analytic expressions can be 
readily worked out. 

As a crosscheck, we also employ the Mathematica package FIRE 48[ and the code 



Apart [49| to perform an independent calculation. Thanks to the integration-by-part (IBP) 
algorithm built in FIRE, it turns out that all the required master integrals (Mis) become 



just the conventional 2-point and 3-point scalar integrals as given in [13j. The final results 
generated by this more automatic approach exactly coincides with those obtained from the 
partial-fractioning method. 

As a third consistency check, the calculation is redone by exchanging the order between 
helicity projection and loop integration. That is, we first utilize the programs Apart and 
FIRE at the amplitude level, which are more cumbersome and time-consuming, yet still 
technically feasible. Once the IR-finite T-matrices are obtained, we then project out each 
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intended helicity amplitudes at the very end. We again find the exact agreement with the 
previous two methods. This calculation can be viewed as a strong support for the validity 
of the 4-dimensional helicity projectors given in Appendix [A] 

Each individual diagram in Fig. [21 being ultraviolet convergent, albeit contains logarith- 
mic infrared divergence. Dimensional regularization (DR) is adopted to regularize those 
IR singularities. Upon summing all the diagrams, the ultimate expression for each helicity 
amplitude becomes IR finite, which endorses the validity of NRQCD factorization approach 
for these exclusive T decay processes. 

Following f lTUj) . we express the reduced helicity amplitude in the three-gluon channel as 



a 



J 

3g;X,X 



(iV c 2 -l)(iV c 2 -4) al 



mi 



57T mi 



Ami 



f 



A,A V 



(14) 



The color factor reflects the fact the three "Abelian" gluons in Fig. [2J must bear an odd C- 
parity, thus proportional to d abc d abc = (N% — l)(A r ( 2 — 4)/iV c , where d abc denotes the totally 
symmetric structure constants of SU(N C ) group. 

All the loop effects are encapsulated in the complex-valued dimensionless functions 
/^-(r). Their full expressions are somewhat lengthy, so will not be reproduced here 5 . 
On the other hand, the profiles of these functions over a wide range of r are explicitly shown 
in Figs. El SI El 



(0,0) 




(1,0) 



10"' 0.25 




10"' 0.25 



FIG. 3: Real and imaginary parts of f®^( r )- The solid curves correspond to the exact results, and 
the dashed curves represent the asymptotic ones taken from (|15a ) and (|15b ). The vertical mark is 
placed at the phenomenologically relevant point r = 0.10. 



In our previous calculation reported in Ref. prior to performing the loop integration, we erroneously 
carried out the Dirac trace in 4 spacetime dimensions. This is an unfortunate mistake which contradicts 
the spirit of DR. In the current work, we take all the Lorentz vectors (both loop and external momenta) 
as the D = 4 — 2e dimensional objects when calculating the Dirac trace. However, we would like to stress 
that, the helicity projectors listed in Appendix El which are derived by simply assuming D = 4, are still 
applicable in this situation. That is because the ultimate amplitudes arc UV, IR finite, so it does not 
matter whether the external momenta are taken as 4 — 2e- or 4-dimensional in the intermediate steps. 
Finally, we note that the analytic expressions of the various reduced helicity amplitudes in the three-gluon 
channel for T — > J/ip + Xo0,i,2 markedly differ from those given in Ref. [34( 1 ■ 
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(0,0) (0,1) (1,0) 




FIG. 5: Real and imaginary parts of /^( r )- The solid curves correspond to the exact results, and 
the dashed curves represent the asymptotic ones taken from (fToT ) to (fT5]i ). 

Theoretically, it is interesting to ascertain the asymptotic behaviors of these reduced 
helicity amplitudes in the limit mj, ^> m c . As mentioned before, we anticipate to see the 
logarithm scaling violation to the naive power-law scaling given in ([8]). The asymptotic 
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expressions of the functions read (J = 0, 1, 2): 
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(15j) 



For readers' convenience, all the asymptotic results of f^{f) are also shown in Figs. El 



IU El in juxtapose with the corresponding exact results. We observe that for most helic- 
ity configurations, the asymptotic results do not coincide well with the exact ones at the 
phenomenologically relevant point r = m 2 /m^ ~ 0.10. 

From f )15p , one confirms that the scaling violation is indeed of the logarithmic form. More 
interestingly, we see that the occurrence of the double-logarithm In 2 r is always affiliated 
with the helicity-suppressed (|A + A| > 0) decay channels. This is compatible with the 



empirical pattern observed in Ref. |32|, |50J. Indeed, such double logarithms have previously 
also been observed in the helicity-suppressed bottomonium exclusive decay processes, e.g. 
T — > J/ijj + Tj c [13], i] b — » J/ipJ/ip [4]. The light-cone approach is presumably the proper 
tool to handle these process-dependent double logarithms. Unfortunately, due to some 
longstanding problems, it remains as a challenge to have a systematic control over these 
double logarithms appearing in the NRQCD short- distance coefficients |50|. 

Finally we mention one peculiarity associated with the helicity channel T — > J/i[>(\ = 
±1) + x c i(A = 0). Recall that at LO in a s , this helicity amplitude from the single-photon 
channel is suppressed by an extra factor of r with respect to the HSR, as can be seen in 
( 113b[) . Nevertheless, from ffl5H ). we find that this helicity amplitude in the 3g channel just 
possesses the correct power-law scaling as dictated by the HSR. This implies that the power 
suppression of the LO single-photon amplitude is purely accidental. 
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FIG. 6: Some representative lowest-order diagrams that contribute to T — > jgg —> J /if) + XcJ- 



C. One-photon-two-gluon channel 

From the ratios of the measured inclusive T decay rates in three different channels, as 
listed in fl2J), an educated guess is that the exclusive decay channel T — > jgg — > J /if) + \ c j 
yields the least important contribution. Nevertheless, for the sake of completeness, let us 
finally assess the contribution from this radiative decay channel. 

Similar to the strong decay channel, this radiative decay channel also starts at the one- 
loop order. Some typical LO diagrams have been illustrated in Fig. [6j For simplicity, and to 
be commensurate with the approximation adopted for the single-photon channel in Sec. lIII Al 
we only retain those photon-fragmentation diagrams, whereas the neglected diagrams are 
the identical as Fig. [2] except with the gluon outside the loop replaced by the photon. We 
wish that these fragmentation-type diagrams constitute the dominant contributions, which 
is certainly the case for the transversely-polarized J/if). 

Following the steps outlined in Sec. IIII B[ one then projects out the 10 required helicity 
amplitudes, with all the polarization vectors (tensors) of T, J/if> and Xcj eliminated from 
the loop integral. Nevertheless, it appears to be much more difficult than in the three-gluon 
channel to employ the partial fraction to simplify the encountered one-loop integrals. 



Fortunately, the powerful Mathematica packages FIRE |48j and Apart [49| can still 
be successfully applied to reduce the general higher-point tensor one-loop integrals into a set 
of Mis. With the aid of the IBP algorithm built in FIRE, all the involved Mis become just 
the standard 2-point and 3-point scalar integrals. The analytic expressions of these scalar 



integrals can be found in [13|, |5l|, whose correctness have been numerically verified by using 



LoopTools 47 



Analogous to the three-gluon decay channel, each individual diagram in Fig. [H] is UV 
finite but IR divergent. After summing all the diagrams, the ultimate expression for each 
helicity amplitude turns out to be completely IR finite. 

In accordance with f fl0|) . the reduced helicity amplitude in the radiative decay channel 
can be expressed as 



a 



J 



1 ei,e c aag 



Alt 



mi 



Ami 



(16) 



The inclusion of an extra factor r~' A ' is reminiscent of the photon fragmentation mechanism: 
when the J/if> becomes transversely polarized (A = ±1), the corresponding helicity amplitude 
would receive a 1/r enhancement with respect to the nominal HSR. 

All the loop effects are encoded in the complex- valued, dimensionless functions g~(r(r). 
Their full expressions are somewhat lengthy, and will not be reproduced here. On the other 
hand, the profiles of these functions over a wide range of r are shown in Figs. [S] and 
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FIG. 7: Real and imaginary parts of g® ?(r). The solid curves correspond to the exact results, and 

A. A 

the dashed curves represent the asymptotic ones taken from (|17a ) and (|17b ). The vertical mark is 
placed at the phenomenologically relevant point r = 0.10. 
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FIG. 8: Real and imaginary parts of g\ Ar). The solid curves correspond to the exact results, and 

A, A 

the dashed curves represent the asymptotic ones taken from (fTTh ) to (fT7e ). 

It is curious to know the asymptotic behavior of the reduced helicity amplitudes in this de- 
cay channel. After some manipulations, we find the asymptotic expressions of the functions 
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FIG. 9: Real and imaginary parts of g^ ^{r). The solid curves correspond to the exact results, and 
the dashed curves represent the asymptotic ones taken from (fTTl O to (jTTj ). 
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g J- (J = 0, 1,2) to read 
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In Figs, d El |9l we also juxtapose these asymptotic results of gi\{f) with the exact 
results. For most helicity configurations, the asymptotic results seem not to converge well 
with the exact ones for the phenomenologically relevant point r = m 2 /m 2 ~ 0.10. 

A quick survey on (|17|) reveals that that the scaling violation is again of the logarithmic 
form. More interestingly, the same pattern of double logarithms still holds: the occurrence 
of the double-logarithm In 2 r is always affiliated with the helicity-suppressed decay channels. 

We close this section by making a simple observation. As can be seen from Fig. |8l the 
imaginary parts of the functions g 1 1- are generally nonzero, though some of which vanish 
asymptotically. At first sight, this may contradict Landau- Yang theorem because Xci 99 
should be strictly forbidden. This implies that, though the apparent two-gluon cuts in Fig. [6] 
do not contribute, the other cuts that simultaneously pass through the bottom and charm 
quark lines must yield nonvanishing contributions to the imaginary parts for T — > J/tp + Xci- 

IV. PHENOMENOLOGY 



We are now in a position to make concrete predictions for the decay rates of T — > 
J/*P + Xco,i,2, by plugging flT2}, flU and (PJ into (03). 

In the numerical analysis, we take the various bottomonia and charmonia masses from the 



2012 PDG compilation [l[: M r 



"(15) 



9.460 GeV, M T 



(25) 



10.023 GeV, M T 



(35) 



10.355 
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GeV, Mj/^ = 3.097 GeV, M Xc0 = 3.415 GeV, M Xcl = 3.511 GeV, M Xc2 = 3.556 GeV. 
These inputs are used to determine |P| according to (J4]), which appears in the phase space 
factor in ( TTTj) . To calculate the squared matrix elements, particularly the functions /^~(r) 

and 9 J x j( r )i we instead adopt the following values for the quark masses: m c = 1.5 GeV, 
rrib = 4.7 GeV, corresponding to r = 0.102. Later we also wish to predict the branching 
fractions for the decays T(nS) — >■ J ' /ip + Xco,L2 ( n — 1; 2, 3). For this purpose, we take the 
total decay width of various T states from '% T[Y{IS)} = 54.02 keV, r[T(25)] = 31.98 
keV, and T[T(3S)] = 20.32 keV, respectively. 

For simplicity, we fix the electromagnetic fine structure constant as a = 1/137. For 
the strong coupling constant, any value evaluated with the renormalization scale ranging 
from 2mb to 2m c seems to be acceptable. This scale ambiguity constitutes one of the most 
important sources of theoretical uncertainty. Without much prejudice, we simply choose a 
medium value a s (mb) = 0.22. 

As for the wave functions at the origin for various quarkonia, we use the values obtained 
from the Buchmuller-Tye potential model Q: | J R T(15) (0)| 2 = 6.477 GeV 3 , | J R T (2S)(0)| 2 = 
3.234 GeV 3 , \R T{3 s){0)\ 2 = 2.474 GeV 3 , |i?j/*(0)| 2 = 0.81 GeV 3 , |^ c(Ui2 (0)| 2 = 0.075 GeV 5 . 



TABLE I: The values of the reduced helicity amplitudes a 7 , a 3g and a igg associated with each 
helicity channel for T(IS') — > J/tp + Xc0,i,2- For simplicity, a factor of 10 -4 has been pulled out of 
each entry. 
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In Table HI we tabulate the values of the reduced helicity amplitudes a 7 , a 3g and a igg 
associated with each helicity channel for T(IS') — > J/ij) + Xco,i,2- O ne sees that the relative 
phases among each amplitude vary channel by channel, and particularly there seems no 
universal phase pattern between the single-photon and three-gluon channel. Our finding 



contradicts the universal relative phase conjecture made in [53[ . 

One curious question is whether the relative strength between three distinct decay mech- 
anisms bears the roughly same pattern as that in T inclusive decay, as given in ([2]). A 
very crude guess is that each helicity amplitude for this exclusive decay process might be 
proportional to \/B mc i, so that one naively expects 

|a 7 | : |a3 9 | : |a 795 | ~ 1 : 3.3 : 0.5, (18) 

for each helicity configuration. 
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TABLE II: The polarized and the polarization-summed partial widths and the corresponding 
branching fractions, rlfi , and rl^ g represent the individual decay rates from the single-photon, 
three-gluon, and one-photon-two-gluon channels affiliated with the process Y(1S) — > J/ip + Xcj 
(J = 0,1,2), respectively, while rj^{ denotes the decay rate which incorporates all three decay 
mechanisms according to (|lip . The two rightmost columns also give the unpolarized decay rates 
and branching fractions by summing over all possible helicity configurations. All the partial decay 
widths are in units of eV. 
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Inspecting Table [H we find that for many helicity configurations, the three distinct helicity 
amplitudes do exhibit the similar hierarchy as given in ( JT8l) . i.e., \a,3 g \ > \a 7 \ > \a igg \, 
though the radiative decay amplitude are often much more suppressed. However, there 
are also a few notable exceptions, e.g., for the polarized decays T — > J/ip(±l) + Xci(0) 6 , 
T — > J/tp(±l) + x c2 (±l), T — > J/-0(±1) + x c2 (±2), the magnitude of the radiative decay 
amplitudes is comparable, or even greater, than that of the single-photon amplitudes; for 
the decay T — > J/ip(0) + Xc2(0), the single-photon amplitude is even greater in magnitude 
than the respective three-gluon amplitude. 

In Table [TTl we also list the polarized decay widths for Y(IS') — > J/ip(X) + Xco,i,2(A) 
for each independent helicity configurations (A, A), together with the polarization-summed 
results. To readily visualize the interference effect among three distinct decay mechanisms, 
we also tabulate the individual decay rates from the single-photon, three-gluon, and one- 
photon-two-gluon channels, respectively, as well as the full decay rate given by (TTTj) . Table [TTl 
reveals that a 7 and in T(IS') — > J /ip + Xco,2 are subject to destructive interference. 

The hierarchy among the polarized decay widths in different helicity configurations, from 
either an individual decay mechanism or the complete contributions, seems hardly to obey 
the HSR as indicated in ([8]). The pattern is abnormal for T — > J/V + Xco, where the helicity- 
suppressed (±1, 0) channel even possesses a bigger decay rate than the helicity-favored (0, 0) 



6 The unnaturally small single-photon amplitude in this channel is due to the accidental suppression factor 
received by c} , as can be seen in (I13bl) . 
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TABLE III: The unpolarized partial decay widths and the corresponding branching fractions that 
incorporate all three distinct decay mechanisms. The superscript [J] characterizes the correspond- 
ing decay process T(15, 25, 35) — >• J/ip + Xcj (J = 0, 1, 2). 
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channel; for T — > J/ip+Xc2, the largest polarized decay rates are associated with the helicity- 
suppressed configurations (±1,0) and (±1,±1), whereas the smallest polarized decay rate 
is associated with the HSR-favored (0, 0) state 7 . These symptoms can presumably be 
attributed to the fact that the mass ratio m c jm^ rs 1/3 might not be small enough to 
warrant the asymptotic counting rule. 

Finally in Table IHIl we tabulate our predictions of the partial decay widths, which are 
calculated according to ([7j), together with the corresponding branching fractions, for the 
various processes T(nS) — >■ J/ip + Xco,i,2 (n = 1,2,3). All the three decay mechanisms 
are incorporated. We observe that the decay branching fractions satisfy the ordering 8 
ffl 1 ] > g[o] ^ g[2]^ w it n the first two reaching the order of 10 -6 . We note that all these 
predicted branching ratios are compatible with the various experimental bounds on J/ip or 
Xcj inclusive production rates in T(1S,2S) decays, as given in (TjQ). 

As a simple consequence of the LO NRQCD prediction, there should exist a 77% rule in 
the hadronic decay of the T system, in analogy with the famous 12% rule in the ip system: 

B[T(2S) -» hadrons] _ B[T(2S) -> e+ej _ ( , 

B[T(1S) -> hadrons] B[T(1S) ->• e+e ) K } 

Not surprisingly, the ratios of B[T{2S) -»■ J/^ + Xcj] to fi[T(lS) -»■ J/ip + Xcj] in Table HH 
are indeed compatible with this rule. 

Thus far, the Belle experiment has collected about 102 million T(1S) samples and 158 
million T(25) samples. According to Table [LTT| the Belle experiment is expected to have 
produced about 130 and 170 T{1S,2S) -> J/ip + Xco events, 500 and 650 T{1S,2S) ->■ 
J/ip + Xci events, 20 and 30 T(1S, 2S) — >■ J/ip + Xc2 events, respectively. 

Experimentally, there are two possible methods to detect the J/ip + Xd0,i,2 signals. The 
first is to reconstruct both the J/ip and Xc0,i,2 events. The clean and copious decay modes 
of Xa,2 are the El radiative transitions Xa,2 J/ip ± 7, with the branching fractions of 
34.4% and 19.5%, respectively jlj. The J/ip meson can be most cleanly tagged through 
the leptonic decays into e + e~ and with combined branching ratios about 12% |l(. 



7 Note this situation is quite different from the continuum production process e + e~ — > J/ip(X) + x C 2(A), 
where the (0,0) and (±1,0) channels make the dominant contributions to the unpolarized production 
cross section [32j j . 

8 It is interesting to compare the relative importance of the various exclusive production channels J/ip + \cj 
(J = 0,1,2) in T decay and continuum production. For e + e~ — > J/ip + Xc0,i.2, one finds that the 



production rate for J/ip + Xco is about one order of magnitude greater than those for J/ip + Xd,2 3l[ l32j|. 
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Although this method has the advantage of bearing very low background level, taking into 
account the reconstruction efficiencies for two J /if) and one photon, one may end up with 
too few signal events to be practically useful. 

The second method is to only reconstruct one J /if) — > event, then fit the recoil mass 
spectrum against the J /if) to estimate the number of Xco,i,2 peak events. This method will 
not depend on the concrete decay modes of Xco,i,2- F° r l° w statistics of signal events like in 
our case, this method is much more superior to the preceding one. As a matter of fact, this 
method has already been used by the Belle collaboration to impose the upper bound for 
the exclusive bottomonium decays X&o,i,2 — >■ J/if)J/if>, J/ifnf)' (H3 |. 

In fitting the recoil mass spectrum of J /if), the net detection efficiency for T(1S,2S) — > 
J/i 3 + Xeo.1,2 is estimated to be around 4% (similar for all Xeo.1,2), with the reconstruction 
efficiency for J/if> — > included [55| . Therefore, the numbers of the observed T(1S, 2S) — > 
J/if> + Xci events are expected to be 500 x 4% = 20, and 650 x 4% = 26, respectively. Since 
only one J /if) is reconstructed, the background level in real data may not be very low. With 
only 20 reconstructed signal events, it seems quite challenging for the signal significance to 
reach the 5a level, and a larger data pool is needed in order to draw a definite conclusion. 
In the prospective Super B factory, with a luminosity 50 times greater than the current 
B factory, it seems very promising that the decays T(1S, 2S, 3S) — > J /if) + XcO,i wm be 
eventually observed. 



V. SUMMARY 



In this paper, we carry out a comprehensive investigation on the exclusive J /if) + Xc0,i,2 
production in T decay in the NRQCD factorization framework. We have explicitly consid- 
ered three distinct decay mechanisms, i.e., the strong, electromagnetic and radiative decay 
channels. Although there has not yet appeared a rigorous proof on the validity of NRQCD 
factorization approach to these types of double-charmonium production processes, the ex- 
plicit verification for the cancelation of IR divergences in our calculation is rather supportive 
of the positive answer. Moreover, our explicit calculation further supports the previous claim 
that the double logarithms appearing in the one-loop NRQCD short- distance coefficients are 
always affiliated with the helicity-suppressed channels [32j, [50 . 

The branching fractions for the various polarized and unpolarized decay channels 
T(nS) —7- J /if) + Xeo,i,2 {n = 1,2,3) are predicted by incorporating all three distinct decay 
channels at the lowest order. In our case, the relative phase among these decay channels arise 
from the short-distance loop effect, which can actually be calculated in perturbation theory. 
There appears no universal interference pattern, but the three-gluon and the single-photon 
amplitudes often tend to be destructive. We find that T(nS) decays into J /if) + Xci have 
the largest branching fraction, about a few times 10 -6 ; and the decays T(nS) —±J/if) + Xc2 
have the smallest decay branching ratio, only of order 10~ 7 . The current statistics at Belle 
is on the margin of observing these decay channels. If the prospective high-luminosity 
e + e~ facilities such as the Super B experiment can dedicate more machine time on the first 
three T resonances, it should be an ideal place to discover their exclusive decay modes into 
J/if> + XcO,l- 
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Appendix A: Various helicity projectors for T — > J/ifi + Xc0,i,2 

During this work, we have utilized the various helicity projectors to expedite projecting 
out the corresponding helicity amplitudes associated with T — > J /if + Xco,i,2- This helicity 
projection technique, which has already been applied in our previous work on the 0(ct s ) 



correction to the process e + e~ — > J /if) + Xc0,i,2 |32|, can significantly reduce the amount of 
labors required for the loop diagram calculations. In this Appendix, we collect the explicit 
formulas for the 10 helicity projectors used in this work. 

Firstly, it is convenient to introduce the transverse metric tensor: 



P^Pu 
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IPI 2 
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P»Pu 
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|P| 2 


M\ 


p|2 



9fiv+ /l/f2 I'D 1 2 (Pf-Qv QfiPv) A/ j- 2 |p|2 ' 



where Q, P and P stand for the four-momenta of T, J/if> and Xcj, respectively. This 
symmetric tensor satisfies the transversity condition g± /JlV P ,J ' = g^^P^ = 0. It further has 
the properties = 2, g ±lia gl u = g±^ a g au = g^- 

The decay amplitude A for the process T — > J /if) + Xco can be parameterized as .A = 
A° v e T ej^,(A), where ex and ej/$ denote the polarization vectors of T and J /if), respectively. 

The helicities of J /if) and Xco are labeled by A and A (trivially A = 0). There are only 
two independent helicity amplitudes for this process, which can be deduced by acting the 
corresponding helicity projectors on the amputated amplitude: A° x ^ = W^-A^, up to an 

immaterial phase. The two helicity projection tensors for T — > J /if) + Xco are 

o,o | P |2^ M T ^^J\Mj /f Mr M r J 
n,o = ~9±,u, (A2b) 

where gj_p U is defined in ( IA1I) . These two projectors are normalized as IPo.O;/^ o = 1> 
Pi.o^fPio = 2> anc ^ orthogonal to each other: P 0j o ;At j/ Pi o = 0. 

The decay amplitude A 1 for the decay process T — > J /if) + Xci can be expressed as 
A 1 = Aj lua e T ey / ,^(\)e™ l (\), where A and e Xcl denote the helicity of the xa state and the 
corresponding polarization vector. There are three independent helicity amplitudes for this 
process, which can be deduced by acting the corresponding helicity projectors on the ampu- 
tated amplitude: A l x =• = P^* A^ ua . The three helicity projection tensors for T — > J /ip + Xci 
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read: 
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which are subject to the normalization conditions F i;flua Fj Ua = with i, j signifying one 
of the three helicity configurations. 

Similarly, for the decay process T — > J/ip + Xc2, we can identify the amputated amplitude 
through A 2 = A 2 ll/a i 3 € T €*j l/ ,^(X)e^ (A), where A and e Xc2 represent the helicity of the Xc2 
state and the corresponding polarization tensor. There are in total five independent helicity 
amplitudes, which can be obtained by acting the corresponding helicity projectors upon the 



amputated amplitude: A 2 r 
for T ->■ J/ip + Xc2 as 



P' 



A, A 



We construct the five helicity projection tensors 



1,0 



r 0,l 



v^|P|2 \Mj/^M T v M r 



P -±*Q 

M 2 



9±aB + 



X 



M Xc2 M T 
1 



2y/G 9±l " 
1 / Q-P 



My 



<7±a/3 + 



M Xc2 M T 
Q-P 



My 



(A4a) 



pp Vm Xc2 m t a 



P -^Q 



My 



M Xc2 M T 



P, 



M 



2V/2IPI 2 \Mj^M r My 



9± 



Q-P 

M Xc2 M r 



M 



Xc2 



M r 



My 



(A4b) 



" M| V " 



M Xc2 M T 



^ (<7± ^ g± a/3 — g± na g± up — g± ^ g±ua)- 



(A4c) 

, (A4d) 
(A4e) 



These five projection operators satisfy the normalization conditions F^ap f^ a/3 = ^Sij (i, 
j corresponding to one of the five helicity configurations), except that Po,o ; 



^0,0 



I. 



All the helicity projectors are derived by obeying the exact decay kinematics. Never- 
theless, in this work we only target at the LO accuracy in v expansion. Therefore, Upon 
applying these projectors to infer the intended helicity amplitudes, it is eligible to make 
the following substitutions for the various quarkonium masses: M-^^s) ~ % m b { n — 1)2,3), 



and M 



M 



Xc0,l,2 



2m c . Implementing these approximations considerably simplifies the 
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corresponding loop calculation. 
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